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1.  INTRODUCTION 

This  work  is  a theoretical  study  of  tlie  equilibrium  properties'"^ 
of  a relativistic  electron  beam.  The  major  recent  interest  in  the 
equilibrium  properties  of  relativistic  electron  beams  originates  from 
several  diverse  research  areas.  These  include;  (a)  fundamental 
experimental  and  theoretical  studies  of  equilibrium  and  stability 
properties  of  magnetically  confined  nonneutral  plasmas  in  mirror, 
cusped,*'  and  uniform'^  magnetic  field  geometries,  (b)  research  on 
collective  - effect  accelerators  that  utilize  the  intense  self  fields 


1.  Davidson,  R.  C.,  Theorv  of  Nonneutral  Plasma  (Benjamin,  Reading, 

MA. , 1974) . * 

2.  Bennett,  W.  H.,  "Magnetically  Self-Focussing  Streams,"  Phys.  Rev. 
45,  1934,  p 890. 

3.  Hammer,  D.  A.,  and  Rostoker,  N.,  "Propagation  of  High  Current 

Relativistic  Electron  Beams,"  Phys.  Fluids  1970,  p.  1831. 

4.  Gratreau,  P.,  "Generalized  Bennett  Equilibria  and  Particle  Orbit 
Analy’sis  of  Plasma  Columns  Carrying  Ultra-High  Currents,"  Phys. 
Fluids  2J.,  1978,  p 1302  . 

5.  Davidson,  R.  C.,  and  Uhm,  H.  S.,  "Thermal  Equilibrium  Properties  of 
an  Intense  Relativistic  Electron  Beam,"phvs.  Fluids,  in  press  1979. 

6.  Armstrong,  C.  M.,  Hammer,  D.  A.,  and  Trivelpiece,  A.  W.,  "Non- 
neutral Electron  Plasma  V'lith  a 2.5  MeV  Mean  Energy,"  Phys.  Lett. 
A^,  1976  , p 413. 

7.  Davidson,  R.  C.,  and  Uhm,  H.  S.,  "Influence  of  Strong  Self- 

Electric  Fields  on  Ion  Resonance  Instability  in  a Nonneutral 
Plasma  Column,"  Phys.  Fluids,  1977,  p 1938. 

8.  Davidson,  R.  C.,  and  Uhm,  H.  S.,  "Influence  of  Finite  Ion  Latmor 

Radius  Effects  on  the  Ion  Resonance  Instability  in  a Nonneutral 
Plasma  Column,"  Phys.  Fluids,  1978,  p 60. 

9.  Striffler,  C.  D.,  Kapetnnakos,  C.  A.,  and  Davidson,  R.  C.,  "Equili- 
brium Properties  of  a Rotating  Nonneutral  E Layer  in  a Cusped 
Magnetic  Field,"  Phys.  Fluids,  18,  1975,  p 1374. 


4 


NSWC/'<JOL  TR  78-182 


of  an  electron  cluster, the  electron  cyclotron  resonance  mode  (the 
Auto  Resonant  Accelerator) , ' ' and  the  converging  waveguide  , ' (c)  re- 

search on  the  high-power  microwave  generation’’”'^  and  (d)  studies 
of  electron  beam  propagation  through  a neutral  gas  or  bacl^ground 
plasma. ° Although  these  research  areas  have  different  goals  and 
objectives,  they  have  in  common  th€  need  to  understand  the 


10.  Destler,  W.  W. , Uhm,  H.  S.,  Kim,  H.,  and  Reiser,  M.  P.,  "Study 
of  Collective  Ion  Acceleration  in  Vacuum"  J.  Appl.  Phys.  in 
press  (1979 ) . 

11.  Sloan,  M.  W.  , and  Drummond,  '-7.  E.,  "Autoresonant  Accelerator 
Concept,"  Phys.  Rev.  Lett.  3^,  1973,  p 1234. 

12.  Sprangle,  P.,  Drobot,  A.,  and  Manheimer,  W. , "Collective  Ion 
Acceleration  in  a Converging  Wave  Guide,"  Phys.  Rev. Lett.  36 , 

1976,  p 1180. 

13.  Granatstein,  V.  L.,  Sprangle,  P.,  Parker,  R.  K.,  Pasour,  J., 
Herndon,  M. , Schlesinger,  S.  P,,  and  Setter,  J.  L.,  "Realization 
of  Relativistic  Mirror;  Electromagnetic  Backscattering  from  the 
Front  of  a Magnetized  Relativistic  Electron  Beam,"  Phys.  Rev. 

AU,  1976  , p 1194. 

14.  Flyagin,  V.  A.,  Gaponov,  A.  V.,  Petelin,  M.  I.,  and  Yulpatov, 

V.  K.,  "The  Gyroton,"  IEEE  Trans.  Microwave  Theory  and  Techniques, 
MTT-25,  1977,  p 514. 

15.  Uhm,  H.  S.,  and  Davidson,  R.  C.,  "Quasi-Linear  Theory  of  Cyclo- 
tron Maser  Instability,"  Phys.  Fluids  22,  in  press  (1979)  and 
references  therein. 

16.  Reiser,  M.  "Laminar-Flow  Equilibria  and  Limiting  Currents  in  Magnetically 
Focused  Relativistic  Beams",  Phys.  Fluids  477  (1977). 

17.  Kligman,  R.  L.,  "Electric  Conductivity  in  a Beam,  Plasma  System," 
NSWC/WOL  TR  77-146,  15  Sep  1977. 

18.  Briggs,  R.  J.,  Hester,  R.  E.,  Lauer,  E.  J.,  Lee,  E.  P.,  and 

Spoerle^n,  R.  L.,  "Radial  Expansion  of  Self-Focused,  Relativistic 
Electron  Beams,"  Phys.  Fluids,  1976,  p 1007. 

19.  Lee,  E.  P. , "Kinetic  Theory  of  a Relativistic  Beam,"  Phys.  Fluids 
1^,  1976,  p 60. 

20.  Lampe,  M.,  "Energy  Moment  Method  for  Straggling  of  an  Ultra- 
Relativistic  Electron  Beam,"  NRL  Memorandum  Report  3221,  Feb 
1976. 
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equilibrium  properties  of  a relativistic  electron  beam  that  is  char- 
acterized by  intense  self  electric  and  self  magnetic  fields. 

The  theoretical  analysis  of  relativistic  electron  beams  is  carried 
out,  neglecting  the  atomic  process  and  discrete  particle  interactions 
(i.e.,  binary  collisions),  because  the  collective  processes  are 
assumed  to  dominate  on  the  time  and  length  scales  of  interest.  The 
equilibrium  properties  are  examined  within  the  framework  of  the 
Vlasov-Maxwel 1 equations.’  The  analysis  is  carried  out  for  an  infi- 
nitely long  electron  beam  as  shown  in  Figure  1.  The  mean  radius  of 
the  electron  beam  is  denoted  by  a.  The  mean  motion  of  the  electron 
i>eam  is  in  the  azimuthal  and  axial  directions,  producing  the  self 
magnetic  fields.  Cylindrical  polar  coordinates  (r,  0,  z)  are  in- 
troduced with  the  z axis  along  the  axis  of  symmetry;  r is  the  radial 
distance  from  the  z axis, 

and  0 is  the  polar  angle  in  the  plane  perpendicular  to  the  z axis. 

The  relativistic  electron  beam  is  immersed  in  a cylindrical  back- 
ground plasma  and  is  partially  charge  neutralized  by  extra  positive 
ions.  We  also  expect  that  the  background  plasma  has  an  additional 
equal  number  of  positive  ions  and  electrons.  The  positive  ions  are 
assumed  to  be  immobile  (m^  . Thus  we  have 


n^(r)  = ri.  (r)  + n^(r)  , 
e b p 


(1.1) 


0 , . 0 , , , 0 , , 
ni(r)  = ' 

where  n?(r)  and  n^(r)  are  the  equilibrium  beam  and  plasma  densities, 
0 ^ 0 

respectively,  n^Cr)  and  n^(r)  are  the  total  electron  and  ion  densi- 
ties, respectively,  and  f^=const  is  the  fractional  charge  neutrali- 
zation. It  is  also  assumed  that  the  space  current  density  generated 
by  the  beam  electrons  is  partially  neutralized  by  the  return  current 

of  plasma  electrons,  with  the  fractional  current  neutralization  f = 

m 

const. 

Equilibrium  properties  of  the  electron  beam  are  calculated  for 
the  class  of  rigid-rotor  equilibria  described  by  kinetic  equilibrium 
distributions  of  the  form  f^(H,  P^,  P^)  = f^(H  - a)^Pp,P^),  where  H 
is  the  total  energy,  P^  is  the  canonical  angular  momentum,  P^  is  the 
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axial  canonical  momentum,  and  is  the  canonical  anqular  velocity. 

It  is  worthwhile  to  emphasize  that  in  rigid-rotor  equilibrium,  the 
electron  distribution  function  depends  exclusively  on  the  linear 
combination  of  the  energy  H and  the  canonical  angular  momentum  P^. 
One  of  the  characteristic  features  of  the  rigid-rotor  equilibrium  is 
that  the  mean  azimuthal  velocity  a beam  electron  fluid  element  is 
rigid-rotor,  i.e.,  ( r ) =u)j^r  , where  V^j^Cr)  is  the  mean  azimuthal 

velocity  of  a beam  electron  fluid  element.  Equilibrium  properties 
are  investigated  for  a variety  of  choices  of  f^(H-u)j^Pp,  P^)  to  de- 
termine the  sensitivity  of  density,  temperature,  and  self-field 
profiles  to  beam  injection  conditions. 

To  make  the  theoretical  analysis  tractable,  the  following  addi- 
tional simplifying  assumptions  are  made  in  describing  the  electron 
beam  equilibrium  by  the  steady-state  (9/3t=0)  Vlasov-Maxwell 
equations ; 

(a)  The  equilibrium  properties  are  independent  of  z (9/3z=0)  and 
aximuthally  symmetric  (9/90=0)  about  the  z axis. 

(b)  It  is  further  assumed  that  the  motion  of  the  beam  electrons 
is  predominately  in  the  axial  direction,  and  that  the  transverse 
momentum  is  small  in  comparison  with  the  characteristic  axial 
momentum , i.e., 

Pj.  -*■  Pf,  ''<9^  > (1.2) 

where  p = (p^.  Ppr  P^^  Is  the  particle  momentum  in  cylindrical  co- 
ordinates. Consistent  witli  Fq.  (1.2),  we  simply  omit  the:  axial  self 
magnetic  field  in  the  subsequent  analysis,  since  the  influence 

of  the  axial  self  magnetic  field,  produced  by  the  mean  rotation 
Vg(r)  = Wj^r  of  the  electron  beam,  on  the  equilibrium  properties  is 
negligibly  small  in  comparison  with  the  influence. of  azimuthal  self 
magnetic  field  B^(r).  A simple  estimate  gives 


,/“b»  \ 

V “b=  \ V ? ) 

zb  0 \ zb  / 


<<  1, 


(1.3) 
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where  is  the  mean  axial  velocity  of  a beam  electron  fluid 

element . 


Jentral  to  a Vlasov  description  of  electron  beam  equilibria  are 
the  single-particle  constants  of  the  motion  in  the  equilibrium  fields 


E°(y)  = E°(r) 


e = - 
~r 


3r 


?r 


(1.4) 


and 


B°(x)  = Bg  4 Bg{r) 

3A®(r) 

= "o  !z  - 3?^ !e 

where  e^,  ^0,  and  are  unit  vectors  in  the  radial,  azimuthal,  and 
axial  directions,  B^  is  a uniform  applied  magnetic  field,  is 

the  equilibrium  electrostatic  potential  and  A^{r)  is  the  z-  com- 
ponent of  equilibrium  vector  potential.  For  azimuthally  symmetric 
equilibrium  with  3/3z=0,  there  are  three  single  particle  constants 
of  the  motion.  These  are  the  total  energy  H, 

„ ,24^22,1/2  - , ,,  ,, 

H=(mc  4cp)'  - e;Q(r),  (1.6) 

the  canonical  angular  momentum  P„ 

Pe  = ^(Pe  - j 

and  the  axial  canonical  momentum 


P_  = p - ^^(r) 
z c z 


(1.8) 


where  c is  the  speed  of  light  in  vacuo , -e  and  m are  charge  and  rest 
mass  of  electron,  respectively,  and  = eB^/mc  is  the  nonrelati- 

visitic  electron  cyclotron  frequency  in  the  uniform  applied  magnetic 
field  Bq . Note  that  without  the  applied  magnetic  field  (8^=0),  the 
term  proportional  to  in  Eq.  (1.7)  vanishes. 

The  radial  electric  field  E^(r)  and  azimuthal  self  magnetic  field 
Ba(r)  are  determined  self-consistently  from  the  Maxwell  equations 
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and 


13  3 ^ 


(r) 


47re  ( 1-fg)  (r)  , 


(1.9) 


1 i_  3_ 

r 9r^9r 


A^(r) 


4iTe--  _ . 0,  .,,0  , > 


(1.10) 


where  and  are  the  fractional  space  charge  and  current  neutral- 
izations respectively,  of  the  system  [see  Eq.  (1.1)].  Here,  n^(r) 

0 ' ° 
and  V , (r)  are  the  beam  density  and  axial  velocity  profiles  defined 

in  terms  of  f,  (H-oj.Pq,?  ) by 
b b o z 

n”(r)  - /d^p  . (1.11) 

n»(r)V»b(r)  = /d^p  £» (H-c^P^ ,P^) , (1.12) 

where  = p^c'^/(m'^c  +c'^p'^) . Equations  (1.6)  - (1.12)  can  be 
used  to  investigate  a wide  variety  of  beam  equilibrium  properties 
in  the  following  two  sections.  Moreover,  the  mean  azimuthal  velocity 
Vgb^’^^  can  also  be  obtainecj  from 

2 2 4 2 '>  1/2 

where  Vg  = pgC /(m  c + c p")  ' 


The  organization  of  this  report  is  the  following.  In  Section  2, 
we  investigate  the  equilibrium  properties  of  rigid-rotor  electron 
beam  with  uniform  axial  drift  velocity.  The  equilibrium  distribution 
function  for  rigid-rotor  electron  beam  with  uniform  axial  drift 
velocity  is  given  by 


fb  - f^"-“b''6-eb'=''z> ' 


which  depends  on  H , Pg,  and  exclusively  through  the  linear  com- 
bination H-a)j^Pg-6j^cP^ . Here  BjjC  is  a constant.  The  equilibrium 
properties  of  a stiff  electron  beam,  which  is  characterized  by  zero 
axial  temperature  and  is  described  by  the  distribution  function  of 


NSWC/WOL  TR  78-182 


the  form 

is  investigated  in  Section  3.  Here,  G(H-<.^p  ) is  an  arbitrary 
function  of  and  is  a constant. 
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2.  ELECTRON  BEAM  EQUILIBRIA  WITH  UNIFORM  AXIAL  DRIFT 

We  first  consider  the  class  of  rigid-rotor  Vlasov  equilibria  with 
uniform  axial  drifts  where  f^  depends  on  H,  and  exclusively 
through  the  linear  combination  H-u)|^P^-Bj^cP^, where  =1-1/7^  =const. 
That  is,  we  consider  beam  equilibrium  distributions  of  the  form 

'b  = • '2.1) 

We  will  find  that  equilibria  described  by  Eq.  (2.1)  correspond  to 
rigid-rotors  with  angular  velocities  and  with  uniform  axial 
velocities  ( r)  =Bj^c=const . In  this  regard,  it  is  useful  to 
transform  variables  in  Eq.  (2.1)  from  momentum  variables  (p^.,  Pg , P^) 
appropriate  to  the  laboratory  frame  to  momentum  variables  (p^,  Pg , p^ 
appropriate  to  a frame  of  reference  moving  with  velocity  Bj^c  e^ 

The  relevant  transformation  is  given  by 


Pr=Pr'Po"P0'Pz"^b^Pz'^''  ’ + BjjD^/mc)  , (2.2) 


where  >=  ( 1+P^/m 'c^)  and  > ' = ( 1+p^/m^c^) 


After  some  straight- 


2 2 2 2 2 2 
forward  algebra  with  the  approximation  p'  =(p^  ♦•Pg  c , we 


find 


^e-'b‘"^^"cb/2) 


b b 


(2.3) 


where  to  . =a)  , /y.  =eB„/Y.  me  is  the  characteristic  electron  cyclotron 
cb  cb  b O b 

frequency  and  i|^(r)  is  defined  by 


2.  _2 


i()(r)=(l/2>Vj^m(u)j^-j0^jj-u)j^)  r -e  ( >^Q-B^3A^) 

2 2 

= Cl/2)Yj^m(u)j^u)^j^-u)j^)  r -e4'(r). 


(2.4) 
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and  't’(r)  is  the  effective  self-field  potential 


(2.5) 


In  order  to  evaluate  the  macroscopic  equilibrium  properties 
defined  in  Eqs . (1.11)  - (1.13),  it  is  also  necessary  to  consider  the 

transformation  of  d^p  to  beam  frame  variables  d^p'.  It  is  straight- 
forward to  show  from  Eq . (2.2)  that  d^p  transforms  according  to 

P ' 

(2.6) 


> ' me 


Substituting  Eqs.  (2.1)  - (2.6)  into  Eq . (1.11),  we  find  that  the 

equilibrium  beam  density  is  given  by 


n^(r)=4TTYb/dp"p"^f® 


...2 


iYbm 


+ >p  (r) 


(2.7) 


where  the  integration  variable  p"  is  defined  by 


P"^  = p;^(p'-Y^,m(.^^r)^p;2. 


Without  loss  of  generality,  we  assume  that  the  effective  self-field 
potential  4'(r)  vanishes  at  r^O,  i.e.,  't'(r=0)=0.  Similarly  the 


function  ^(r)  also  vanishes  at  r=0.  Defining  nj^=n^(r=0),  and 


F(ii^)  = 


® 2 2 
/dp"p"^f°  (n^£-  + — + h 


/dp"p”^f^(^.  ^) 


(2.8) 


2Ybm 


Eq.  (2.7)  can  be  expressed  as 


0 


nb(r)  = n^  F[iJ;(r)  ] . 


(2.9) 


It  is  also  straightforward  to  show  from  Eqs.  (1.12),  (2.1)  and 
(2.3)  that 


H 


‘ i 
H 

4 


ii 
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n°(r)vO,(r)  = /d^p  ^ fO 


= 6j^cnj^(r) 


(2.10) 


We  note  from  Eq.  (2.10)  that  uniform  over  the  beam  cross 

section  with 


V^j^(r)  = 0j^c  = const. 


(2.11) 


This  is  a characteristic  property  of  electron  beam  equilibria 
described  by  an  electron  distribution  function  which  depends  on  H 
and  exclusively  through  the  linear  combination  H-gj^cP^.  More- 
over, substituting  Eq.  (2.3)  into  Eqs . (1.13)  and  (2.1),  we  readily 
find 


(2.12) 


It  is  obvious  from  Eq.  (2.12)  that  the  mean  aziumthal  velocity  of 
the  beam  electron  fluid  element  is  that  of  a rigid-rotor.  Again, 
this  is  a typical  characteristic  of  the  rigid-rotor  beam  equilibria 
in  which  the  distribution  function  also  depends  on  H and  P^ 
exclusively  through  the  linear  combination  of  H-Wj^Pg.  We  therefore 
conclude  from  Eqs.  (2.11)  and  (2.12)  that  any  electron  beam  equili- 
brium described  by  a distribution  which  is  exclusively  a function  of 
H-a)bp0“Sb^^z' ^ ^ a rigid-rotor  equilibrium  distribution  with  a uniform 
axial  drift. 

The  remaining  important  macroscopic  quantity  that  characterizes 
the  electron  beam  is  the  equilibrium  temperature  profile  of  the  beam 
electrons.  The  effective  beam  temperature  T^(r)  is  determined  from 


n°(r)T°(r)4/d^p[VrPr+(v0-V°b)  (Pe-<Pe^) 


^ ^'^z-'^zb’  (Pz-<Pz>)  1 fb  V^z^ 


(2.13) 


where  v^  =p (m^c^+c^p^)  , and  <Pj>=  ( /d^pp^  f^)  / ( /d^pf^) 
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Substituting  Eq.  (2.3)  into  Eqs.  (2.1)  and  (2.13)  gives 


n“(r)T»(r) 


roc 


(2.14) 


2 2 2 

where  use  has  been  made  of  the  approximation  p'  <<m  c . One  of  the 
important  features  in  Eq.  (2.14)  is  that  the  effective  beam  tempera- 
ture T^(r)  is  isotropic,  i..e.,  ( r ) ( r)  =T®j^  ( r)  =T®  ( r)  , where 

"^rb^^^'  etfective  radial,  axiumthal,  and 

axial  beam  temperatures,  respectively,  defined  by 


(Vj-<v.  >)  (p.-<Pj>)  f°  (2.15) 

where  (P f ^ 1 / ( /d^pf^l  . Examples  of  anisotropic 

beam  temperatures  will  be  discussed  in  the  next  section. 

The  Maxwell  equation  for  the  effective  self-field  potential 
4’(r)  can  be  obtained  from  the  linear  combination  of  Eqs.  (1.9)  and 
(1.10) , and  we  have 


I If  p Ir 


(2.16) 


2 "'2 

where  aipj^=4TTnj^e  /Y^^ni  is  the  characteristic  plasma  frequency-squared 
and  use  has  been  made  of  Eqs.  (2.9)  and  (2.11).  It  is  evident  from 
Eqs.  (1.9),  (1.10)  and  (2.5)  that  the  potentials  4)g(r)  and  A®(r)  are 

related  to 


1-f. 


♦o(r)  = 


1-f  -6;:+f  8. 
e b mb 


rrFTT  ' 


(2.17) 


and 


A®(r)  = V^— T 

l-VV^m^b 


(2.18) 


Therefore,  once  the  effective  potential  4'(r)  is  determined  from 
Eq.  (2.16)  the  equilibrium  self  fields,  E^  (r)  =-9((iQ/9r  and 
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s s 

Bg (r) =-3A^/9r , can  be  calculated  self-consistently  from  Eqs.  (2.17) 
and  (2.18). 

Eliminating  the  effective  potential  't'(r)  in  Eq . (2.16)  in  favor 
of  the  function  v(r)  and  defining  a new  dimensionless  function 


ip(r)  = ip(r)/T  , 


(2/19) 


where  T = const  is  a measure  of  the  thermal  energy  spread  for  the 
distribution  function  f^,  we  find  the  nonlinear  differential  equation 
for  I)/  ( r)  : 

13  3 , , , 2, 

r 37  3r*"^'  “ — <“cb"b-“b' 


”T 


E^(l-f_-Bd+f  6nF[i^(r)  ] 


^e  "b 


(2.20) 


The  measure  of  thermal  energy  spread  T is  a parameter  characteristic 
of  the  beam  temperature.  Once  F is  known,  the  equilibrium  problem 
has  been  reduced  to  solving  the  nonlinear  differential  equation  (2.20) 
for  i^(r).  It  is  useful  to  express  Eq.  (2.20)  as 


— r i)/  ( r ) = - 

r 7r  ^ ' T 


2Ybm 


+ a)^a)'[l-F(i;/)  ] I , (2.21) 

+ 

where  the  laminar  cold-fluid  rotation  frequencies  (*1^  are  defined  by 

“b'  ' ir“=b*'"cb-''>^^^]  - <2-22> 


in  which 


K = 4u)^u)~  = 2u)^.(l-f-6^+fB^) 
b b pb  e b mb 


(2.23) 


is  a measure  of  the  combined  strengths  of  the  equilibrium  self- 
electric field  (1-f^)  and  the  equilibrium  self-magnetic  field 
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[ 


[ 

i 

t 

I 


From  the  definition  of  ij)  and  a careful  examination  of  Eq.  (2.21), 
we  find  that  (r=0  ) =0=  OC'/'lr)  ^_q  . For  physically  acceptable  solu- 
tions of  Eq , (2.21)  corresponding  to  radially  confined  beam 

equilibria  (with  n^(r)  approaching  zero  for  sufficiently  large 
radius  r),  F (. ) is  a monotonic  decreasing  function  of  increasing 
\p  with  F(t1'=0)=l  at  i(/(r=0)=0,  by  Eq . (2.8),  and  ijj(r)  is  a monotonic 

increasing  function  of  r.  An  algebriaic  manipulation  shows  in  a 
straightforward  way  that 


1 

r Tr 


r 


3 

JF 


-I 

V I 

r=0 


^^b"’  + 


= (2.24) 

We  note  from  Eq.  (2.24)  that  for  a radially  confined  equilibrium, 
the  beam  rotation  frequency  is  restricted  to  the  range 

-'b'  1 % 1 ' ^2.2L>; 

which  is  only  a sufficient  condition  for  radial  confinement.  The 
complete  necessary  and  sufficient  condition  for  radially  confined 
equilibrium  will  be  obtained  in  the  specific  example.  For  in  the 
range  of  Eq.  (2.25),  we  find  [r~^  ( 3/3r)  (r3i(//3r)  ) . Therefore, 

beginning  with  ij;(0)=0,  if  we  integrate  Eq . (2.21)  for  small  r and  , 

we  find  that  4j(r)is  monoton  i ca  1 1 y increasing  near  the  origin.  However, 
for  large  r,  ij;(r)  is  not  always  to  be  a monotonic  increasing  func- 
tion of  r.  Hence,  the  inequality  in  Eq . (2.25)  is  a sufficient  bi.t 

not  necessary  condition  for  a radially  confined  equilibria. 

In  Eq.  (2.21),  is  a variable  (as  yet  unspecified)  rotation 
parameter.  Moreover,  the  nature  of  the  solution  for  t|;(r)  exhibits 
a strong  dependence  on  the  value  of  relative  to  the  values  of 
and  , which  are  the  laminar  cold-fluid  rotation  frequencies 
for  a beam  with  constant  density  fij^.  For  future  reference,  we  plot 
versus  K [Eq.  2.23)]  in  Figure  2.  The  segment  in  the  right-half 


i- 

i 


I 
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plane  in  Figure  2 corresponds  to 

K > 0),  (2.26) 

i.e.,  the  repulsive  space-charge  force  on  a beam  fluid  element  exceeds 
the  magnetic  focussing  force  associated  with  the  azimuthal  self 
magnetic  field.  On  the  other  hand,  the  left-half  plane  in  Figure  2 
corresponds  to 

6^(l-f  ) > 1-f  (or  K < 0) , (2.27) 

b m e 

in  which  case  the  self-magnetic  force  on  a beam  fluid  element  exceeds 

the  repulsive  space-charge  force.  We  note  from  Figure  2 that  the 

u)£  curves  are  defined  in  the  entire  left-half  plane,  and  converge  to 

the  common  point  co."  = uj  /2  in  the  right-half  plane  for  a value  of 

beam  density  satisfying  -B^+6.  f ) = o)  , . We  also  note  from 

^ ^ pb  e b b m cb 

Figure  2,  and  Eqs.  (2.22)  and  (2.24)  that  no  equilibrium  solutions 
exist  in  the  region  where 

The  reason  is  simply  that  equilibrium  space  charge  force  is  so  strong 
in  this  region  that  the  total  magnetic  attractive  forces  (self-field 
plus  applied  field)  cannot  overcome  the  repulsive  electrostatic  force. 
Moreover,  it  is  emphasized  that  without  the  applied  magnetic  field, 

Eq.  (2.22)  can  be  further  simplified  as 

2u)^  = (-K)^^^  , (2.2n 

and  the  co^  curves  are  defined  only  in  the  left-half  plane  of 
Figure  2. 

To  complete  the  general  analysis  of  this  equilibrium  configuration, 

we  define  the  effective  beam  radius  a by  the  relation 

00 

Aj^a^  =2/  drrn^(r).  (2.30) 

2 2 

Moreover,  we  introduce  Budker's  parameter  v defined  by  v=Nj^(e  /me  ), 
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where 


= 2n/  drrn°(r)  = a' 


(2.31) 


is  the  number  of  electrons  per  unit  axial  length  of  the  beam.  For 
a general  profile  n®(r),  we  find 


2 2 
u)  ,3 
pb 


(2.32) 


which  is  one  of  the  most  important  parameters  in  discussing  electron 
beam  equilibriun..  Evidently,  v/Yj^  ^ 1 accordingly  as  a < 2c/u)pj^. 

We  now  analyze  Eq.  (2.20)  for  several  specific  choices  of  distribution 
function  f^,  which  illustrate  specifically  the  qualitative  behavior 
described  above. 

2 .A  Sharp  - Boundary  Beam  Equilibria 


As  a first  example,  we  consider  a choice  of  distribution  function 
for  which  the  beam  density  is  identically  zero  beyond  some  radius 
r^,  i.e.,  n®(r>r^)=0.  In  particular,  we  consider  the  distribution 

2 

U(T+^  -H+o).  P^+6,  CP^) 
n ^b  b 9 b z 

fO  = (2.33) 

-„t„bP9+6bCP^>^^2 

b 

where  the  Heaviside  step  function  U(x)  is  defined  by 


U(x)  = 


1 , X i:  0 , 


0 , X < 0 


(2.34) 


The  distribution  in  Eq.  (2.33)  is  chosen  in  this  subsection  simply 
because  the  nonlinear  differential  equation  of  il'  in  Eq.  (2.20)  becomes 
linear  and  can  be  solved  analytically  for  the  entire  range  of 
interesting  physical  parameters.  Substituting  Eq.  (2.33)  into  Eqs. 
(2.8)  and  (2.9),  and  carrying  out  the  straightforward  algebra,  we 


find 


F(4»)  = (l-i|/(r)]  U (l-H^) 


(2.35) 
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and  the  electron  density  profile 


l-!^(r) 


n.  (r)  = h. 


for  ij'(r)<l, 
for  (r)  >1 . 


(2.36) 


Before  determining  the  solutions  of  the  nonlinear  differential 

equation  of  Eq . (2.20),  it  is  convenient  to  treat  separately  the 

cases  (a)  1-f  >6^(1-!  ),  and  (b)  8^(l-f  ) >l-f„. 

e b m b m e 


(a ) Equilibrium  Dominated  by  Space-Charge  Forces 
2 2 

(1-f  ^3w-f  B ).  We  first  consider  the  region  of  the  right- 
e b mm 

half  plane  in  Figure  2,  where 

0 i K = < '‘^cb'  (2.37) 

In  order  to  reduce  the  number  of  independent  parameters  in  Eq . (2.20), 

we  introduce  the  dimensionless  variable 


X = r/d^  , (2.38) 

where  d^  is  defined  by 

dg  = 16T/>j^mK  . (2.39) 

Substituting  Eq.  (2.35)  into  Eq . (2.20),  and  making  use  of  Eqs . (2.38) 

and  (2.39)  gives 


where 


1 3 


8; 


/ 8 (1-ij;)  , 0<r<r  , 

^ 0 , r>r^  , 


(2.40) 


(2.41) 


and  r^  is  the  confinement  radius  satisfying  ilj(r^)=l.  For  t<r^,  the 
solution  to  Eq.  (2.40)  is 


4^(X) 


(;-l)  [I„(2/TX)-1]  , 


(2.42) 
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where  Iq (x)  is  the  modified  Bessel  function  of  the  first  kind  of 
order  zero.  In  obtaining  Eq . (2.42),  use  has  been  made  of  il/(r=0)=0 

[see  the  discussion  following  Eq.  (2.23)). 

Since  the  function  li'(r)  must  be  greater  than  zero  for  all  r,  it 
is  quite  obvious  from  Eq . (2.42)  that 

C > 1 (2.43) 

is  the  necessary  and  sufficient  condition  in  the  case  for  a radially 
confined  beam  equilibrium.  [Note  that  Iq  (“)  = “. 1 , Making  use  of 
the  definition  of  r.  in  Eq.  (2.41),  Eq . (2.43)  can  be  expressed  as 

Eq.  (2.25)  for  a positive  value  K.  Thus,  we  conclude  from  Eq,  (2,43) 
that  the  space-charge  dominated  electron  beam  described  by  Eq . (2.33) 

is  in  radially  confined  equilibria  provided  the  beam  rotation  fre- 
quency is  in  the  range  u)“< . Substituting  Eq,  (2.42)  into 
Eq.  (2.36),  we  find  the  electron  density  profile 


;-(;-!)  Iq(2,TX) , 0<r<r^  , 


n-(r)  = 


r>r 


where  the  beam  confinement  radius  r is  obtained  from 

c 


Iq  (2r^r^/dQ)  = ;/(;-!). 


(2.44) 


(2.45) 


The  effective  beam  radius  a defined  in  Eq,  (2.30)  can  be  simply 
expressed  as 

a^  = - (^7/2)  (C-l)r^dQl^(2/7r^./dQ)  (2.46) 

In  Eq,  (2.46)  Ij^(x)  is  the  modified  Bessel  function  of  order  one. 

Note  from  Eqs.  (2.45)  and  (2.46)  that  the  beam  confinement  radius  r^ 
and  the  effective  beam  radius  a increase  to  infinity  when  ; approaches 
unity.  Moreover,  the  measure  of  the  thermal  energy  spread  T plays  a 
significant  role  in  determining  the  confinement  radius  r^  and  beam 
radius  a [see  Eqs.  (2.39),  (2.45)  and  (2.46)). 
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(b)  Equilibrium  Dominated  by  Se 1 f-Maqnetic  Forces  (6.  -f  6.  >l-f 

~ “ ■■  ■ D in  O € 

We  now  consider  the  left-half  plane  in  Figure  2,  where  ' " 


K — 4 uj,  ^ 0 • 

b b 


(2.47) 


2,2 

In  this  region,  B.-f  6,  > l-f-*  i.e.,  the  self-magnetic  force  on  a 


b mb  "e ' 

beam  fluid  element  exceeds  the  repulsive  space-charge  force.  n,  is 
convenient  to  analyze  Eq.  (2.20)  with  the  radial  variable  r normal- 
ized in  units  of  the  Bennett  pinch  radius  b^  [see  Eq . (2.83)].  De- 
fining 


X"  = r/bp. 


(2.48) 


where 


b^  = - d^  = 16T/>.m;K|, 


(2.49) 


Eq.  (2.20)  can  be  expressed  in  the  form 


1 3 


/ 8 (l-i/<)  , 0 vr<r^  , 


X"  - -8^;  + 


(2.50) 


0 , r<r  , 

' c ' 


where  ; is  again  defined  in  Eq.  (2.41).  For  the  solution  v(X") 


that  satisfies  iij(X"=0)=0  is  given  by 


«|;(X")  = (l-0  [1-Jq  (2vTX")  ] , 


(2.5i: 


where  dQ(5c)  is  the  Bessel  function  of  the  first  kind  of  order  zero. 
Once  again,  in  order  to  have  a positive  ;|;(r)  for  all  r,  it  is 


necessary  that 

C < 1 , 


(2.52) 


which  also  can  be  expressed,  as  in  Eq.  (2.25),  as  u),  for  a 

D o b 

negative  K value.  In  this  case,  the  inequality  in  Eq.  (2.25)  is 
only  a sufficient  condition  for  radially  confined  equilibrium. 

However,  Eq.  (2.25)  is  not  a necessary  condition  for  radial  con- 
finement, since  the  function  i(/(r)  in  Eq.  (2.51)  is  always  smaller 


1 


c 
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where  the  right-hand  side  inequality  is  obtained  from  C { ^/2 v 2 ) > 1 , 
are  also  very  interesting.  In  terms  of  the  beam  rotation  frequency 
Eq.  (2.54)  can  be  expressed  as 


than  unity  for  a certain  range  of  the  physical  parameter  l,  even  if  ; 
satisfies  Eq . (2.52).  In  order  to  have  radially  confined  equilibrium, 

the  value  of  the  function  \l  (r)  must  be  greato'-  than  unity  for  r r^ , 
where  r^  satisfies  v(f^)=l.  In  finding  a necessary  condition  for 
radial  confinement,  it  is  important  to  investigate  the  properties  of 
Eq.  (2.50)  and  the  Bessel  function  ^^(x)  which  has  the  minimum  value 


-0.402  at  .x=a. 


Here  aQ2^=3.83  is  the  first  root  of 


Jl(aoi)=0  Jj^(x)  is  the  Bessel  function  of  order  one. 


Since  the 


value  of  i)i(r)  is  always  greater  than  unity  for  r r^  when  i>0  [Examine 
carefully  Eq . (2.50)1,  we  find  the  necessary  and  sufficient  condition 


(2.53) 


for  a radially  confined  equilibrium. 

The  equilibrium  properties  of  the  electron  beam  corresponding 
to 

0 - 1 1 - (2.54) 


(2.55) 


where  is  determined  by 


In  this  case,  the  value  of  the  function  ;J^(r)  would  be  less  than  unity 
for  several  subranges  of  r.  Thus,  the  density  profile  would  have 
subsidiary  peaks  at  successive  radii  and  the  central  beam  would  be 
surrounded  by  coaxial  cylindrical  beams  for  in  the  range 
0 - ■ t ( '■'q ])  1 physical  significance  of  such  a configura- 

tion is  not  yet  clear.  We  shall  call  this  condition  tlie  multi  annular 
beam  equilibrium. 

Substituting  Et] . (2.51)  into  Eq . (2.30),  wc  find  tlie  density 
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profile  of  beam  electrons 

n°(r)  = A^^[c+(1-c)Jq  (2/7r/bQ)  ] , 


(2.56) 


for  either  radially  confined  or  oscillatory  equilibrium.  For  a 

radially  confined  equilibrium,  in  which  the  beam  rotation  frequency 

0)^  satisfies  Eq.  (2.53),  the  beam  density  becomes  zero  at  the 

radial  confinement  radius  r . On  the  other  hand,  if  cj,  is  in  the 

c b 

ranye 


< CO,  < 

b b b 


or  CO, 


<^b  ""^b' 


(2.57) 


we  find  from  Eq.  (2.51)  that  ip(r)  [and  hence  n^(r)]  exhibits 
oscillatory  behavior  with  its  maximum  value  always  less  than  unity. 
Evidently  from  Eq.  (2.56),  in  this  oscillatory  equilibrium,  the 
beam  density  n^(r)]  decreases  from  at  r=0 , to  as  r-^<». 

Moreover,  the  beam  density  exhibits  an  oscillatory  dependence  on  r, 
with  characteristic  radial  wave  number 


(2.58) 


where  b^  is  the  Bennett  pinch  radius. 

The  effective  beam  radius  a is  estimated  by  substitution  of 
Eq.  (2.56)  into  Eq . (2.30)  and  we  find 

a^  = !;r^+(/I/2)  (l-c)r^bQj^(2/Tr^/bQ)  , (2.59) 

where  the  ratiial  confinement  raiiius  r is  obtained  from 

c 

jQ(2/2r^/bQ)  = ^/(C-1).  (2.60) 

When  the  self-magnetic  force  exceeds  the  repulsive  space-charge 

force,  the  maximum  radial  confinement  radius  r of  the  electron 

cm 

beam  described  by  Eq.  (2.33)  is  given  by 

^cm  " “01^0^^*^^  ' (2.61) 
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which  corresponds  to  :,  = 1 - ( 1- Jq  ( ^ ) 1 ~ ^ . Moreover,  the  radial 
confinement  radius  r^  (hence  the  effective  beam  radius  a)  increases 
by  increasimj  the  mesaure  of  thermal  energy  spread  T [see  Eqs.  (2.49), 
(2.59)  and  (2.60) ) . 

A summary  of  the  equilibrium  properties  for  the  electron  beam 

described  by  Eq . (2.33)  is  presented  in  Figure  3.  Figure  j is 

similar  to  Figure  2 except  for  the  addition  of  a plot  of  <2^  versus 

K = in  the  left-half  plane  where  the  self-magnetic  force  exceeds 

the  replusive  space-charge  force.  Wo  reiterate  that  for  an  electron 

->  o 

beam  dominated  by  space-charge  force  ( ^~fc-‘'b~  ^m^'b^  ' necessary 

and  sufficient  condition  for  a radially  confined  equilibrium  is  given 
by  Eq.  (2.25).  On  the  other  hand,  Eq . (2.53)  is  the  necessary  and 

sufficient  condition  for  a radially  confined  beam  equilibrium  domi- 
nated  by  a se  I f-maqnpt  i c force(8^-f.^.''j‘^  1-f^^J  . These  are  illustrated 
by  the  cross-hatched  region  in  Figure  3. 

Finally,  we  calculate  the  beam  temperature  profile  defined  in 
Eq.  (2.14).  Substituting  Eq . (2.33)  into  Eq . (2.7)  and  (2.14),  it 

is  straightforward  to  show  that 

T°(r)  = iTll-ii.(r)  ) , (2.62) 

which  decreases  from  T/2  at  r=0,  to  zero  at  for  a radially 

confined  equilibrium.  Evidently,  the  measure  of  thermal  energy  spread 
T is  directly  related  to  the  beam  temperature.  We  therefore  conclude 
from  Eqs.  (2.46),  (2.49)  and  (2.59)  that  the  effective  beam  radius  a 

increases  considerably  when  the  beam  temperature  is  increased. 

2.0  Diffuse  Beam  Equilibrium  - The  Generalized  Bennett  Pinch 

As  a second  example,  we  consider  the  choice  of  distribution 

function 

^b^"'  ^O'^z’  = i exp|-(Il-ui^^P^,-h^^cr^)/T]  , (2.63) 

b 

where  the  normalization  constant  is  defined  by 
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= TT  7^ THTI  expdnc^/Y^^T).  (2.64) 

'b  (2TT>,mT) 
b 

The  equilibrium  properties  of  the  electron  beam  described  by  Eq.(2.63) 

5 

have  been  investigated  in  the  previous  literature  . However,  we 

5 

repeat  some  of  the  arguments  in  the  previous  study  in  order  to  make 
this  report  more  complete.  Moverover,  the  distribution  function  in 
Eq.  (2.63)  is  used  more  often  in  describing  the  electron  beam  pro- 
perties than  any  other  distribution  function,  mostly  because  the 
Gibbs  distribution  in  Eq.  (2.63)  corresponds  to  the  unique  thermal 
equilibrium  distribution  to  which  the  beam  electrons  would  evolve 
through  self  collisions.  Equation  (2.63)  corresponds  to  significant 

2 

generalization  of  the  distribution  function  first  analyzed  by  Bennett 

in  that  Eq.  (2.63)  allows  for  a mean  azimuthal  rotation  of  the  beam 

electrons  when  (oi^/O).  We  will  therefore  find,  for  appropriate  chcice 

of  u)j^,  that  radially  confined  beam  equilibria  exist  even  in  the 

absence  of  a partially  neutralizing  positive  ion  background  that 

weakens  the  repulsive  space  charge  force.  The  reason  is  simply  that 

the  u),  rd  xB„d  force  on  an  electron  fluid  element  is  in  the  inward 
b 0-z 

direction,  and  therefore  acts  as  a confining  force  that  can  overcome 
repulsive  electrostatic  forces. 

From  Eqs.  (2.8),  (2.9),  (2.19),  and  (2.63),  we  find 


F (lij)  = exp  [-4'  (r)  ] 

and  the  density  profile  of  the  beam  electron 
n^  (r)  = expl-t  (r)  1 
The  nonlinear  equation  of  i'{r)  in  Eq . (2. 

13  3 , , ■''b”,  2,  'b"  „ 

F ' T^<“cb"b-“b'  - rr 

where  use  has  been  made  of  the  definition 


20)  can  be  expressed 

exp[-iKr)  1 , 

of  K in  Eq.  (2.23) . 


(2.65) 


(2.66) 

as 


(2.67) 
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( a ) Nece s s ary  and  Sufficient  Conditions  for  Radially  Confined 
Equilibrium 

We  now  make  use  of  Eqs,  (2,21)  and  (2.65)  for  iKr),  to- 
gether with  the  definitions  in  Eqs,  (2.22)  and  (2.23),  to  determine 
the  necessary  and  sufficient  conditions  for  radially  confined 
equilibria,  i,e,,  equilibria  characterized  by 


lim 

If -*-oo 


0 , ^ 
n^(r) 


lim 

If  -kOO 


tij^expl-4'(r)  ] 


0 


or  equivalently 


lim 

If  >co 


i(^(r) 


In  this  regard,  we  investigate  properties  of  the  solution  to  Eq. 

(2.21),  subject  to  the  boundary  conditions  il;(r=0)=0  and  i()(r+“)=“. 

Evidently,  if  4<(r)  is  a monotonic  increasing  function  of  r,  then 

n^  (r)  =hj^exp  (i())  is  a monotonic  decreasing  function  of  r,  decreasing 

from  A,  at  r=0  to  zero  as  r*-"-’.  Since  the  sufficient  condition  for 
b 

radial  confinement  has  been  already  obtained  in  Eq . (2.25)  for 

general  rigid-rotor  equilibria,  we  investigate  here  the  necessary 
condition  for  radial  confinement  particularly  for  the  Gibbs  dis- 
tribution function  in  Eq.  (2.63). 

2 2 

Equilibrium  Dominated  by  Space-Charge  Forces  (l-f^^8^^-8^^f^)  . 

We  first  consider  the  region  of  the  right-half  plane  in  Figure  4 
where  K>0.  To  show  that  Eq . (2.25)  is  also  a necessary  condition 
for  a radially  confined  equilibrium,  let  us  assume  that  is  outside 
the  range  in  Eq.  (2.25)  (either  or  u)j^<a)j^)  and  that  (r-^®)  =“, 

which  corresponds  to  nP(r*®)=0.  This  will  lead  to  a contradiction. 

+ - ° -1 
For  or  we  conclude  from  Eq.  (2.24)  that  [r  (3/9r) 

( r 9ij;/9r ) ] 0 » so  that  iKr)  decreases  to  negative  values  near  the 

origin.  If  i(i(r)  is  to  eventually  assume  a positive  value  so  that 

lij  (r-*-®)  =®,  then  ((((r)  must  pass  through  a minimum  (Figure  5)  at  r=r,  , 

2 2 . ^ 

say,  where  ii'(r^)-,0,  (9>l^/3r)^_^  =0  and  (9  ii'/9r  ) j,_j.  ^0.  Evaluating 
Eq.  (2.21)  at  r=rj^,  we  find  ^ ^ 
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+ - 


[ 1-exp  (-ijj)  ] 


(2.68) 


which  is  manifestly  negative  for  ti'(rj^)<0  and  satisfying 
u)b<oj^.  [Keep  in  mind  K=4.i)^u)”20  ] . We  therefore  conclude  that 
^b^‘^b  contradictory  to  the  assumption  of  a radially 

confined  equilibrium  with  (r-^*)  =+“.  Hence,  a restriction  of  the 
rotational  parameter  to  the  range  specified  by  Eq.  (2.25)  is  both 
necessary  and  sufficient  for  a radially  confined  equilibrium  dominated 
by  space  charge  force.  This  is  illustrated  by  the  cross-hatched 
region  in  the  right-half  plane  of  Figure  4. 

2 2 

Equilibrium  Dominated  by  Self-Magnetic  Forces  (6.-f  8.  >l-f  ). 

— *■ — i 2 b m b e 


We  now  consider  the  left-half  plane  in  Figure  4 where 
K=4(o^u)j^<0 . A determination  of  the  necessary  and  sufficient  con- 
ditions for  radially  confined  equilibrium  solutions  when  K<0  (Figure 
4)  proceeds  in  a similar  manner  to  the  analysis  in  the  case  K>0.  We 

therefore  summarize  the  essential  results.  For  K=4u)!|'a)~<0 , a restriction 

b b 

of  the  rotational  parameter  to  the  range 


0 


< liJ,  < U)  , 

- b-  cb 


(2.69) 


is  both  necessary  and  sufficient  for  the  solution  to  Eq.  (2.67)  to 
correspond  to  a radially  confined  equilibrium  characterized  by 
n^(r-*-“>)=0  and  ^ =+<^\  This  range  is  illustrated  by  the  cross- 

hatched  region  in  the  left-half  plane  of  Figure  4.  Moreover,  for 
K<0  and  0<u)^<l,  the  beam  density  profile  n^(r)  is  a m^  notonic  de- 
creasing funetion  of  r.  Finally,  for  K<0,  if  is  in  the  range 


u)  . <co,  <u)l  or  (ij,  <to,  <0, 

cb  b b b b ' 

we  find  from  Eq.  (2.67)  that  i|j(r)  [and  hence  n^(r)]  exhibits 
oscillatory  behavior. 

To  summarize  the  results,  we  find  that  radially  confined 
equilibrium  solutions  exist  provided  -“<K=4(i)^a)j^<l  and  lies 


(2.70) 
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2 2 

in  the  cross-hatched  reqion  in  Fiqure  4.  Moreover,  for  8.-f  '’1-f  . 

^ _ b m D e ' 

and  in  the  ranqe  or  the  density  profile 

n^(r),  althouqh  a nonincreasing  function  of  r on  the  average,  does 
exhibit  an  oscillatory  radial  dependence.  Furthermore,  we  note  from 
Eq.  (2.67)  that  if  ot  solution  to  Eq.  (2.67) 

is  ii;(r)=0,  which  corresponds  to  uniform  beam  density  for  all  r. 

(b)  Exact  Analytic  Solutions 

Equation  (2.67)  does  permit  exact  analytic  solutions  for 

certain  special  choices  of  the  rotational  parameter  oj.  , or  the 

+ _ ^ 
quantity  that  measures  the  strength  of  the  equilibrium  self 

fields.  Before  determining  these  solutions,  we  find  from  Eq . (2.13) 

that  the  beam  temperature  profile  for  the  Gibbs  distribution  function 

is  given  by  the  isothermal  value 

T°(r)  = T = const.  (2.71) 

We  now  determine  the  exact  solution  to  Eq.  (2.67)  for  two  special 

2 2 

choices  of  equilibrium  parameters:  (a)  1-f  =i.v-f  hi  , which  cor- 

e b m b 

responds  to  a self-field  force-free  equilibrium,  and  (b) 

^b  <'^cb”“’b^  corresponds  to  an  equilibrium  in  which  the 

magnetic  force  associated  with  the  applied  field  and  the  centrifugal 
force  exactly  cancel. 

Self-^ield  Force-Free  Equilibrium.  In  the  special  case  where 


(2.72) 


i.e.,  the  magnetic  and  electric  self-field  force  on  a beam  fluid 
element  exactly  cancel,  Eq.  (2.67)  reduces  to  the  simple  linear 
equation 


1 

r 


3 3 


(r) 


2>bm  2 

"T  ^‘^'cb%“%^  • 


(2.73) 


The  solution  to  Eq.  (2.73)  is 
4'(r)  = r^/a^. 


(2.74) 
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2 (2T/>j^m) 

^0  ix>.  (u)  . -U),  ) 

b cb  b 

Moreover,  the  density  profile  n^(r)=fij^  exp(-i^)  is  given  by  the 
Gaussian  profile 


(2.75) 


Hb  (r)  = exp 


(2.76) 


The  effective  beam  radius  a is  found  by 


(2.77) 


from  Eq.  (2.30)  for  the  self-field  force-tree  equilibrium.  Note  that 
K=0  (Eq.  (2.72)1  corresponds  to  the  vertical  axis  in  Figure  4 and 

is  required  [ Eq . (2.75)]  in  order  for  the  beam  equilibrium 

to  be  radially  confined  with  a^^O. 

Equilibrium  With  ^ second  example  that  has  an 

exact  analytic  solution,  we  consider  the  case  where  the  centrifugal 
and  applied  magnetic  forces  exactly  cancel,  i.e.,  » 01^ 

equivalently 

io.=0  or  u),  =ui  , . (2.78) 

b b cb 

2 2 

From  Figure  4 and  Eq . (2.69),  we  therefore  require  Sj^-f  j^6^>l-f  ^ for 

radially  confined  equilibrium  solutions.  Making  use  of  Eq . (2.78), 

we  find  that  Eq . (2.67)  can  be  expressed  as 


where 


13  ^ I i \ 8 ) i/\ 

F 7F  ^ = 77  I 

^0 


u2  _ 161 

^0  Yj^ml 


8T/Yj^m 

6u-l+f  ) 
pb  b mb  e 


(2.79) 


(2.80) 


The  solution  to  Eq.  (2.79)  is  given  by 


iKr)  = 2«n(l  + r^/bQ)  , 


(2.81) 
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and  the  beam  density  profile  n^(r)  = exp  {-ijd  reduces  to 
0 


_ _ — 

(1+r^/bJ)" 


(2.82) 


which  is  identical  to  that  obtained  by  Bennett  . Substituting 
Eq . (2.82)  into  Eq . (2.30),  we  find  the  Bennett  pinch  radius 


a = b 


(2.83) 


0' 


which  has  been  introduced  into  Eq . (2.48).  The  equilibrium  self 

fields  can  be  calculated  from  Eqs.  (2.4),  (2.17),  (2.18),  (2.19)  and 

(2.81).  We  find 


and 


Ej(r) 


“’pb<l-fe> 


^ 2~r2- 

Ur  /bp 


0(r)  = - 


^ b"'  2 

^ w . i'..  (1-f  ) 
2e  pb  b m 


(2.84) 


(2.85) 


We  note  from  Eq.  (2.80)  that  the  effective  beam  radius  3=bQ  becomes 

large  whenever  i\^(l-f  ) is  closely  tuned  to  1-f  , i.e.,  whenever  the 
b m e 

self-magnetic  pinching  force  almost  cancels  the  repulsive  electro- 
static force. 

( c ) Numerical  Soluti o n s to  Eg.  (2.6  7) 

The  exact  analytic  solution  to  Eq . (2.67)  is  accessible 

only  in  certain  special  limiting  cases  (see  Eqs.  (2.74)  and  (2.81)). 
Here  we  solve  Ec].  (2.67)  numerically  for  a variety  of  system 
parameters  of  experimental  interest.  To  reduce  the  number  of  in- 
dependent parameters,  we  introduce  the  dimensionless  variable 


X'  = r/a^  , (2.86) 

where  a^  is  tlie  effective  beam  radius  defined  for  self-field  force- 
free  equilibrium  (2.75).  Making  use  of  Eq . (2.86),  Eq.  (2.67) 


reduces  to 


Figure  5 


Plot  of  versus  r for 
qnd  radial  confinement. 
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1 

X' 


■>X 


X 


■) 

7x- 


= 4 [ 1-a  exp (- . ) 


(2.87) 


where  a is  defined  by 

K 


n = 


4 U.',  ( ‘A*  I “ 'a’.  ) 
b cb  b 


(2.88) 


Figure  6 shows  a plot  of  the  normalized  ef^^ective  beam  radius- 
squared  defined  by 


= 2 / dX ' X ' exp 


(2.89) 


‘0 


versus  a (Eq.  (2.88)],  which  is  a measure  of  the  beam  density  and 

self-field  strength.  Note  from  Figuie  6 that  the  effective  beam 

radius  a increases  monotonica 1 ly  from  a^  when  a=0 , to  infinity  when 

a = l,  or  equivalently  or  Moreover,  for  a^O  (i.e., 

8r-f  D?  ^1-f  ),  a decreases  from  a=a,,  when  \=0,  to  zero  when  a tends 
b m b e 0 

to  minus  infinity.  Equation  (2.87)  has  been  solved  numerically  for 
various  choices  of  a.  This  is  illustrated  in  Figures  7 and  3,  where 
C'(X')  is  plotted  versus  X'=r/'a^  (Figure  7),  and  the  corresponding 
normalized  density  profile  nj^(X’)  'fij^=exp  [ - . ( X ' ) ] is  plotted  versus 
X'=r/aQ  (figure  8)  for  several  values  of  a.  Wo  emphasize  from  the 
definition  of  a^  in  Eq . (2.75)  that  the  numerical  analysis  of  Eq . 

(2.87)  is  valid  only  when  the  beam  rotational  parameter  is  in  the 
range  ^ • Note  from  Figure  7 that  (i(X')  is  a monotonic  in- 
creasing function  of  X'  for  each  value  of  a.  Moreover,  for  specified 
.X ' , v{X')  decreases  as  a increases,  with  ,(X')=0  in  the  limit  where 
a=l.  We  also  note  from  Figure  8 that  the  curve  for  a=0  corresponds 
to  the  gaussian  density  profile  in  Eq . (2.76). 

In  the  region  where  the  self-magnetic  focussing  force  exceeds 
the  electrostatic  repulsive  force,  i.e.,  K^O,  it  is  convenient  to 
analyze  Eq.  (2.67)  with  radial  v'ariable  r normalized  in  units  of  the 
Bennett  pinch  radius  b^  derived  for  [Eq.  (2.83)]. 

Defining 


X"=r/b 


0 


(2.90) 
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Eq.  (2.67)  can  be  expressed  in  the  equivalent  form 


^ ^ X"  -glTT  iKX")=8  {n+exp[-i(^(X")  1 } , 


(2.91) 


where 

n = - r.  = 4u)j^(u)j^-u)^j^)/K. 


(2.92) 


The  effective  beam  radius  is  calculated  numerically  from 


2 ® 

^=2  /^dX"  X"exp  [-i(^(X")  ] . (2.93) 

2 2 

Figure  9 shows  a plot  of  a /b^  versus  n for  n^O,  i.e.,  for  radially 
confined  equilibria  with  radial  density  profile  decreasing  monotoni- 
cally  as  a function  of  X".  Evidently,  decreases  monotonically 

as  a function  of  n from  the  value  a/bQ=l  when  n=0  (the  familiar 
Bennett  pinch  radius).  Figure  10  shows  a plot  of  (KX")  versus 
X"=r/bQ  obtained  numerically  from  Eq.  (2.91)  for  various  values  of 
n>0  (radially  decreasing  density  profile),  and  n in  the  range  0>n>-l 
(oscillatory  density  profile).  The  corresponding  density  profiles 
are  illustrated  in  Figure  11.  The  special  case  where  ri=0  in  Figure 
11  corresponds  to  the  Bennett  pinch  equilibrium  density  profile  in 
Eq.  (2.82).  We  note  from  Figure  11  that  the  beam  density  profile  is 
an  oscillatory  function  of  X"=r/bQ  for  n in  the  range  -l<n<0,  which 
corresponds  to  in  the  range  w^<u)^<0  or  [see  Figure  4). 


2 . C Monoenergetic  Beam  in  Beam  Frame 

As  a third  example  of  the  uniform  axial  drift  beam,  we 
consider  the  delta-function  distribution 


4 ■ ^ 


(2.94) 


where  N is  the  normalization  constant.  Substituting  Eq.  (2.94)  into 
Eq.  (2.8  ),  and  ma)cing  use  of  Eq.  (2.19),  we  find 


F(ij^)  = (1-iJ;)^''"  Ud-ip) 


(2.95) 
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Figure  11.  Plot  of  Normalized  Beam  Electron  Density  versus  x"  [Eqs.  (2.66)  and  (2.91)] 
for  several  values  of  n. 
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where  U(X)  is  the  Heaviside  step  function  defined  in  Eq.  (2.34). 
The  corresponding  density  profile  of  beam  electrons  is 


n^(r)  = 


( l-()))  , for  \J;<  1 , 

0 , for  ij/^l. 


(2.96) 


Similarly,  the  effective  beam  temperature  T^(r)  is  obtained  as 
T°  (r)  = jTd-i^^)  U(l-ii/)  . 

The  governing  equation  for  ip  is  given  by 


(2.97) 


1 a 9 , 

F ^ - 


^^b"'  , 2, 

T ^‘^cb%"%^ 


^ b"'^  1/2 

(1-4^)  , for  4)<1, 


2T 


(2.98) 


, for  4;>1, 


where  K is  defined  in  Eq.  (2.23).  Since  the  equilibrium  properties 
for  the  delta-function  distribution  is  very  similar  to  the  previous 
two  examples,  we  present  only  the  numerical  solutions  to  Eq.  (2.97), 
skipping  the  possible  analytic  investigation  of  Eq.  (2.97). 

(a)  Equilibtium  Dominated  by  Applied  Magnetic  Forces 
[a)b(u)cb-cOb_^ 

We  first  obtain  the  numerical  solution  for  the  case 


0 <u),  <tk)  , 

- b - cb. 


when  the  applied  magnetic  force  (term  proportional  to  exceeds 

the  centrifugal  force  (term  proportional  to  u)^)  . Introducing  the 
dimensionless  variable  X'  defined  in  Eq.  (2.86),  Eq.  (2.98)  reduces 
to 


(2.99) 


1 9 V • ^ I —A 

9X'  ^ IX'"^ 


4a(l-4/)^/^,4'<l, 


(2.100) 
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Figure  13.  Plot  of  Normalized  Beam  Electron  Density  versus  x'  [Eqs.  (2.96)  and  (2.100)]  for 
Several  Values  of  a. 
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where  defined  in  Eq , (2.88).  Figure  12  is  a plot 

of  ,(X')  versus  X'=r/a  for  several  values  of  a.  The  corresponding 
normalized  density  profile  nj^  (X  ' ) /nj^=  [ 1- v' (X  ' ) 1 is  plotted  versus 
X'=r/aQ  in  Figure  13.  We  note  from  Figure  12  that  the  qualitative 
feature  of  Figure  12  is  very  similar  to  those  of  Figure  7.  We  also 
emphasize  from  Figures  12  and  13  that  the  curves  for  x=Q  correspond 
to  the  self-field  force-free  equilibrium  (K=0).  A plot  of  normalized 
effective  beam  radius-squared  defined  by 

2 V 

~ = 2 ' dX'X*  11-.  (X’ ) 1 ' “ 
a^  0 

versus  . is  presented  in  Figure  14  for  the  3-function  distribution 
function.  Note  from  Figure  14  that  the  qualitative  feature  of 
Figure  14  is  very  similar  to  Figure  6.  We  therefore  conclude  that 
the  equilibrium  properties  of  electron  beam,  whose  rotational 
frequency  is  in  the  range  0'“'b2'‘'cb'  qualitatively  similar  for 

different  distribution  functions. 

2 2 

(b)  Equilibrium  Dominated  Se  1 f-Magnotic  Forces  ( -f  ) 

— i h b e 

In  the  region  where  the  self-magnetic  focusing  force 
exceeds  the  electrostatic  repulsive  force  (i.e.,  K-0),  it  is  con- 
venient again  to  analyze  Eq . (2.98)  with  ratiial  v'ariable  r normalized 

in  units  of  the  Bennett  pinch  radius  b^  defined  in  Eq . (2.83). 

Making  use  of  Eq . (2.90),  Eq , (2.98)  can  be  expressed  as 


1_  3 v..^ 

X"  3X"  3X" 


V(X") 


8n  + 


8(l-i;d 


, li' " 1 , 


(2.101) 


where  n is  defined  in  Eq.  (2.92).  Figure  15  shows  a plot  of  ii'(X") 
versus  X"=r/bQ  obtained  numerically  from  Eq.  (2.101)  for  various 
values  of  n'O  (radially  decreasing  density  profile),  ii  in  the  range 
-0.5<n<0,  and  n in  the  range  -l^n^-0.5  (oscillatory  density  profile). 
The  corresponding  density  profiles  are  illustrated  in  Figure  16.  We 
note  from  Figure  15  that  the  value  of  becomes  exactly  unity  at 

X"=1.4  for  n=-0.5.  Therefore,  we  have  the  central  and  successive 


48 


NSWC/WOL  TR  78-182 


coaxial  annular  beam  densities  for  n in  the  range  -0.5<n>0  [see 

Figure  16  for  n in  the  range  -0.5<n<0and  the  discussion  following 

the  Eq.  (2.55)].  Figure  17  shows  a plot  of  normalized  effective  beam 

radius  - squared  defined  by 
2 

A ^ . OO 


= 2 r 


dX"X"  [l-t|j(X")  ] 


versus  n for  n>0,  i.e.,  for  complete  radial  confinement. 

In  summerizing  this  section,  we  conclude  the  following.  First, 
the  necessary  and  sufficient  condition  for  radial  confinement  is 
given  by  for  K>0  (the  repulsive  electrostatic  forces  exceeds 

the  focussing  self-magnetic  forces)  and  0 K<0  (the  self- 
magnetic  focussing  force  exceeds  the  repulsive  electrostatic  forces) . 
Second,  for  the  sharp  boundary  equilibrium  dominated  by  self- 
magnetic  focussing  forces  (K<0) , we  have  the  multi-annular  beam 
equilibrium  for  the  rotational  parameter  corresponding  to 
0<(;<  [1-Jq  (Oqj^)  for  the  distribution  in  Eq.  (2.33)  and  correspond- 
ing to  0<^<0.5  for  the  6-function  distribution  in  Eq.  (2.94).  Third, 
oscillatory  electron  beam  equilibrium  also  exists  for  a relavent 
rotational  frequency  Fourth,  in  the  rigid-rotor  electron  beam 

with  uniform  axial  drift,  we  also  have  found  that  the  angular 
velocity  as  well  as  the  axial  beam  velocity  is  uniform  over  the  beam 
cross  section  [see  Eqs.  (2.11)  and  (2.12)].  Moreover,  the  effective 
temperature  of  the  electrons  is  isotropic  [see  Eq.  (2.14)].  Finally, 
we  point  out  that  the  analysis  can  be  extended  for  any  other 
equilibrium  distribution  function. 
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. STIFF  ELECTRON  ni’AM 

In  this  section,  wo  invest  icjato  the  equilibrium  properties  of 
the  stiff  electron  beam  that  is  cliaract  or  i zed  by  the  zero  axial 
temperature.  All  the  electrons  have  the  same  axial  canonical 
momentum  so  that  the  il  i s t r i but  i on  function  is 


2 -1  /T 

where  >^=(1-i'q)  *■  is  a constant,  and  G is  an  arbitrary  function 

determined  by  the  injection  condition  of  the  electron  beam.  The  annular 

velocity  of  beam  equilibria  described  by  Fq.  (1.1)  is  uniform  over 

the  beam  cross  section.  However,  the  axial  velocity  V®,  (r)  of  the 

zb 

electron  beam  fluid  element  is  not  uniform  anymore  over  the  cross 

section.  Making  use  of  Kqs.  (1.6)  anil  (1.7)  , and  the  approximation 
2 2 2 2 

p +p  V vm  c , wo  f ind 


H 


u',  r , 
. b 0 


>pmc-  + 


I'i 

I N m 
z 


( r ) T , 


(T.  2) 


where  the  function  x(r)  is  defined  by 

[)^(r)->^lmc"-e>T^(r)  < > .“'fr‘V‘'cb^  2 

x(r)  = i-  , (3.1) 

the  mass  scale  is  defined  by 

>.,(r)mc^  = (m^c'^+c^p^ ) ^ = I m^c"^ + c"  ( ^ -mi',.c+— A^  (r)l  ‘‘1^^  (3.4) 

2.  ^ Z 0 (I  C 7 » 

2 2 2 

Pj^  = p^+ (p^^--^  ^m(o^^r)  Tj^=  const  is  a measure  of  the  transverse 
thermal  energy  spread  and  ii'^j^=en^/mc  is  the  nonrelat  ivi  st  ic  electron 
cyclotron  frequency.  In  obtainirg  Fq . (3.2),  we  have  assumed  that 

A^  ( r = 0 ) =0  = .f^  ( r = 0 ) without  loss  of  generality.  Note  from  Fq  . (1.4) 
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that  Y2(j^=0)=YQr  so  that  x(r=0)=0  in  Eq.  (3.3). 

'Substituting  Eq.  (3.1)  and  (3.2)  into  Eq.  (1.11)  and  defining 


I duG[Y-mc  +U+X (r)Tx 1 

F(X)=  -53 2 ' (3.5) 

/QduG(YQmc  +u) 


we  find  the  density  profile  of  the  beam  electrons  to  be 


"b(^^  = F[x(r)],  (3.6) 

where  is  the  beam  electron  density  at  r=0.  It  is  also  straight- 
forward to  show  from  Eqs.  (1.12),  (3.1)  and  (3.2)  that  the  axial 

mean  momentum  of  a beam  electron  fluid  element  is  approximate^  by 


n^(r)V°b  (r)  = [ 6qC+  (e/Ygmc)  A®  (r)  ] (x)  . 


(3.7) 


The  transverse  temperature  of  the  beam  electrons  is  defined  by 

n°(r)  Tj(r)=  %/d\[v^p^+(Vg-V°j^)  (Pg-<PQ>)]f°,  (3.8) 


where  v-  and  <p^>  are  defined  in  Eq.  (2.13).  Making  use  of  the 
3 3 2 

2 ,2  4,  2^2.%  = Y n'c  , we  can  simplify  Eq.  (3.8) 

approximation  ynic  =(mc+cp)^  'z  ^ j ' 


as 


n^(r)T^(r)  = 27rY2m/Q  duuG  [ Ygnic^+u+ x (r)  Tx ) . (3.9) 

The  axial  temperature  profile  is  obtained  by  substituting  Eq.  (3.1) 
into  Eq.  (2.15)  and  we  find 

T°(r)  = 0 . (3.10) 

Due  to  the  presence  of  the  factor  6 ( P^-Yq^B^c) in  the  equilibrium 
distribution  function,  the  electron  motion  in  the  axial  direction  is 
cold  and  the  effective  axial  temperature  is  zero. 


54 


NSWC/WUL  TR  78-182 


The  equilibrium  properties  of  the  electron  beam  are  determined 
through  Eqs.  (1.9),  (1.10)  and  (3.5)  - (3.7).  We  now  analyze  these 
equations  for  one  specific  choice  of  distribution  fuction  f®,  in  which 
all  electrons  have  the  same  energy  in  the  laboratory  frame.  This  will 
orient  the  reader  and  illustrate  the  equilibrium  properties  of  the  stiff 
electron  beam. 

As  an  example,  we  consider  the  delta-function  distribution 

G(H-WbP0)  = 6 (H-cjj^Pg  - Yomc^-T_^)  , (3.11) 

where  all  the  beam  electrons  have  the  same  energy.  Substituting  Eq. 
(3.11)  into  Eqs.  (3.5)  and  (3.6),  we  find 


F(X) 


for  X<1, 
for  \>1, 


(3.12) 


and  the  electron  density  profile 


. _ ft  ( 


£l| 

0 ( 0, 


for  X(r)il, 
for  X (r) >1. 


(3.13) 


For  notational  convenience  in  the  subsequent  analysis,  we  define  the 
collisionless  s)cin  depth  6 by 


.2  r2/-2 

6 = C /u)  , 

pb 


(3.14) 


? 2 

where  = 4 Trfij^e  /Yq^  is  the  beam  electron  plasma  frequency  squared. 

The  z-component  of  the  vector  potential  A^(r)  is  determined  by  sub- 
stituting Eqs.  (3.  7)  and  (3.12)  into  Eq.  (1.10).  For  yil  we  have 

1-f. 


19  9 y,  S f . 

? 9F^^  - 


(r)  = 4^en^^6o(l-f^). 


(3.15) 


On  the  other  hand,  for  )^''l,  the  vector  potential  A^(r)  is  determined 
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from 


1 3 3-Sj* 

r ?r  “ °- 


(3.16) 


We  assume  that  r is  the  beam  confinement  radius  so  that  X(r  ) = 1, 
c c 


The  solution  to  the  Eq.  (3.15)  and  (3.16),  which  is  continuous  at 


r = r , has  continuous  first  derivative  at  r = r , and  satisfies 


(r=0)  = 0,  is 


,s,  , >0"'^  ^0 

A^(r)  = ^ 


[I^(r  vT:f^/6)-l),  0<r<r^, 


^0 

^0  '"c’ 


II,,  (r„ .0^/6)  + 


rp/lTf 


m 


m 


xCn^  - 1]  ,r>r^,  (3.17) 
c 


where  I , (x)  is  the  modified  Bessel  function  of  the  first  kind  of 


order  ? . The  corresponding  self  azimuthal  magnetic  field  is  given  by 


X^mc^6f,(l-f  )^ 
B|(r)  = -1-^-^^ ^ 


Il(r  /l-fj^/6)  , 0>r>r^, 


(3.18) 


^Il(rc/r^/6),  r>r^. 


Substituting  Eq.  (3.17)  into  Eq.  (3.7)  gives  the  mean  axial  velocity 
profile 


0 

V,b(r)  - 


(3.19) 


where  is  given  by 


>^(r)  = [1+Y^P^  Iq  (r/T^/6)2]'’5 


m 


(3.20) 
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The  electrostatic  potential  (r)  is  determined  from 


(3.21) 


The  exact  solution  to  Eq.  (3.21)  is  not  tractable  analytically; 
however,  the  asymptotic  solution  to  Eq.  (3.21)  can  be  obtained  in 
various  limits  of  practical  interest,  depending  on  the  size  of  the 
self-field  influence. 

3 Tenuous  Beam 

In  circumstances  where  the  beam  density  is  so  low  that 

N e^  1 ^ ^ 

V _ b ^ 1 (•  c 1 ^ , 

__2  T (3.22) 


>n  4 

0 'y^mc 


the  modified  Bessel  function  Iq(x)  can  be  expanded  according  to 
Ip  (X)  = l+x^/4 , giving 


(3.23) 


The  electrostatic  potential  in  Eq.  (3.21)  can  also  be  approximated 
by 


♦o(r)  = 


Substituting  Eqs.  (3.23)  and  (3.24)  into  Eq.  (3.3),  the 
dimensionless  function  \ (r)  is  given  by 

''^0  + - 2 
" 2^ 

where  the  laminar  cold-fluid  rotation  frequencies  are 


(3.24) 


(3.25) 


defined  by 


-b  ' I 4I(»ob''Y„)'-KI^^2,  _ 


(3.26) 


NSWC/WOL  TR  78-182 


in  which 


(3.27) 


The  necessary  and  sufficient  condition  for  a radially 
confined  equilibrium  is  given  by 


‘"b^“b^“b 


(3.28) 


from  Eq.  (3.25).  The  radial  confinement  radius  r^  is  derived 


f rom  X =1  and  given  by 


=r JZi ^1 

LY^m(u)J-u)b)  (<^b““b^ 


1/2 


(3.29) 


Equation  (3.29)  indicates  that  the  beam  radius  r^  increases  when 
the  beam  temperature  T,^  increases.  Moreover,  in  the  event  that 
is  closely  tuned  to  either  or  the  two  laminar  cold~ 
fluid  rotation  frequencies,  the  beam  radius  r^  becomes  very 
large.  It  is  worthwhile  to  note  the  difference  between  the 
radial  confinement  condition  in  the  previous  section  and  that  in 
this  section.  For  the  electron  beam  with  uniform  axial  drift 
described  in  Section  2,  the  necessary  and  sufficient  condition  for 
a radially  confined  equilibrium  dominated  by  the  self-magnetic 
force  (K<0)  is  given  by  0 we  have  oscillatory 
equilibria  for  K<0  and  rotational  frequency  oj.  in  the  range 

a.  ^ 

However,  the  stiff  monoenergetic  electron 


beam  is  radially  confined  even  if  K<0  and  the  rotational 
frequency  is  in  the  range  ‘^<-b*^“b'^‘^b  reason 

for  this  difference  is  simply  that  the  current  source 
n|^  (r)  (r)  in  Eq.  (3.7)  for  the  vector  potential  A^(r)  has  one 

more  additional  term,  which  is  proportional  to  right- 

hand  side  of  Eq.  (3.7).  This  term  ensures  that  the  vector 
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potential  in  Eq.  (3.15)  increases  monotonical ly  by  increasing 

r.  On  the  other  hand,  in  the  electron  beam  with  uniform  axial 

drift,  the  current  source  (2.10)  doesn't  include 

the  term  proportional  to  A®(r),  and  thereby  leads  to  an 

^ + _ 

oscillatory  equilibrium  for  K<0  and  for  “’b  "^“b  ‘ 

3. B Space-Charge  Neutralized  Beam 

As  a second  example  of  finding  an  analytic  expression  of  the 
dimensionless  function  X (r)  in  Eq.  {3.3)for  a monoenergetic  beam, 
we  consider  the  space-charge  neutralized  beam,  i.e.. 


V- 

N 


'e  = '' 


(3. 30) 


which  ensures  the  elimination  of  the  electrostatic  potential  in 
Eq.  (3.3).  Therefore,  in  order  to  have  a radially  confined 
equilibrium,  it  is  required  that  2 

1 “b^"’cb"%)'' 


>z^^^  . 2, 

— > g (r  , 


= 1 + 


2 


(3.31) 


, 1 % 2 
2 c 


for  r>r  , where  r is  the  radial  confinement  radius  satisfying 
c ^ c 

Here  > (r)  is  defined  in  Eq.  (3.4).  We  can 


find  from  Eq.  (3.31)  that  the  function  g(r  ) is  a monotonically 

2 

increasing  or  decreasing  function  of  r accordingly  as 


Since  the  whole  analysis  in  this  paper  is  valid  only  when 

(w.r/c)^  <<1  [see  Eq.  (1.3)],  Eq.  (3.31)  is  also  valid  only  when 
^ 2 

(ai,  r/c)  <<  1 . Note  from  Eq.  (3.31)  that  in  case 
° 2 2 

Wi.  (u)  1.-W,  ) > T,ui,/y  me  , we  always  have  a radially  confined 
b cb  b 1 b'  0 

equilibrium  even  without  the  self-pinching  forces  (i.e.. 
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'' z case  0"'^“ ''u'^  (u'^j^-a'^)  , a radially  con- 

fined equilibrium  can  exist,  provided  there  is  a reasonable  self- 
magnetic  pinching  force.  However,  finding  the  exact  r^  requires  the 
numerical  analysis  of  F.q . (3.31). 

3-C  Equilibrium  with  = 0. 

The  equilibrium  properties  of  electron  beams  with>‘'|^,  = 0 and  described 
by  Fq . (3,11)  have  been  investigated  in  the  previous  literature.^  Tluis, 

we  repeat  here  some  of  the  arguments  made  in  previous  studies.^  In  the 
case  of  u<j^=0,  the  dimensionless  function  on  \(r)  can  be  simplified  to 


\(r)  = 


[ > , (r)  ^ lmc‘'-eC'p  (r) 


(3.32) 


The  density  profile  is  given  bv 


n.  (r)  = [1  + — ]n  r^-r 

b pmc-4T^'  b'  ^c' 


^ ^ ) n,  , r r 

> ^mc  b'  c 


(3,33) 


Substituting  Eq , (3.33)  into  Eq.  (1.9)  gives 


1 

F TF  = 


cvl'.  (r ) 


, r''r_. 


(3. 34) 


The  solution  to  Eq.  (3.34)  is 


^gmc 


|T^(r.  l-f^/.^)-ll  , O^r^r^, 


(3.35) 


c (r^.  1-f^/.'!)  ^n  :^-ll  ,r'-, 

c 


From  Eq.  (3.35),  the  radial  electric  field  E ( r ) =- 3 3r 
can  be  expressed  as 
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r 


I ^ (r  .l-f^/>5)  , 0-r-r^, 


I.  (r  . 1-f  /S) , r-r  , 
r 1 c e c 


(3.^6) 


The  additional  equilibrium  properties  are  obtained  from  Eq . (3.32) 

by  substituting  Eqs.  (3.4),  (3.17)  and  (3.35)  into  Eq . (3.32). 

However,  this  is  not  straightforward  and  requires  some  numerical 
calculation. 

We  conclude  this  section  by  deriving  some  general  properties. 
First  the  perpendicular  temperatuie  of  beam  electrons  defined  in  Eq . 
(3,9)  can  be  expressed  as 


T?(r)  - Tall  - ■.  (r)  ] 


(3.37) 


It  is  also  straightforward  to  calculate  the  total  current  carried 
by  the  electron  beam  for  the  equilibrium  distribution  function 
defined  in  Eq.  (3.11).  Making  use  of  Eqs,  (3.13)  and  (3.19),  we 
can  express  the  magnitude  of  the  total  beam  current  as 


I=27Te/^^  (mc^/e) 


»'T7f 


m 


IT 


‘"'o''  0^1 


1-f 


m 


= 170008 


o'^o 


r yTTT 
c m 


2 8 


amperes 


(3.38) 


when  use  has  been  made  of  mc^/e  = 17000  in  units  of  amperes. 
Evidently,  I in  Eq.  (3.38)  is  sensitive  to  the  value  of 
The  analysis  in  this  section  can  also  be  extended  to  a diffuse 
electron  beam  equilibrium  described  by  the  equilibrium  distribution 
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differential  equations.  However,  an  analysis  of  these  equations 
can  also  be  carried  out  in  a relatively  straightforward  manner. 
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4.  CONCLUSIONS 

In  this  paper,  we  have  investigated  the  equilibrium  properties 
of  a rigid-rotor  relati'.  istic  intense  electron  beam  described  by 
the  equilibrium  distribution  function  of  the  form 
f^(H,P^,r^)  = which  ensures  a rigid-rotor 

equilibrium  with  the  angular  frequency  uij^.  The  theoretical  analysis 
of  relativistic  electron  beams  has  been  carried  out,  neglecting  the 
atomic  process  and  discrete  particle  interactions  (i.e.,  binary 
collisions),  because  the  collective  processes  are  assumed  to 
dominate  on  the  time  and  length  scale  of  interest.  The  positive  ions 
have  been  assumed  to  be  immobile,  providing  partial  space-charge 
neutralization . 

Various  equilibirum  properties  were  calculated  in  Section  2, 
for  a relativistic  electron  beam  with  uniform  axial  drift  velocity. 

A second  order  differential  equation  (2.20)  has  been  obtained  for  the 
dimensionless  function  4'(r)  which  governs  the  equilibrium  properties 
of  the  electron  beam.  We  have  investigated  this  differential 
equation  analytically  and  numerically  for  several  specific  examples 
of  the  electron  distribution  function,  thereby  finding  the  necessary  and 
sufficient  condition  for  a radially  confined  equilibrium.  One  of 
the  important  features  for  rigid-rotor  electron  beam  with  uniform 
axial  drift  is  that  an  oscillatory  equilibrium  exists  for  the 
rotational  frequency  in  the  range  ‘*’b  “^b ^ ' Moreover, 

the  effective  temperature  profile  is  isotropic  for  this  distribution 
function. 

In  Section  3,  the  equilibrium  properties  were  calculated  for  a 
stiff  electron  beam,  in  which  all  the  beam  electrons  have  the  same 
axial  canonical  momentum,  thereby  leading  to  a zero  axial 
temperature.  Particularly,  the  monoenergetic  electron  beam 
properties  have  been  analytically  investigated  in  detail.  The 
expression  for  total  current  carried  by  the  monoenergetic  electron 
beam  was  derived  in  terms  of  the  beam  intensity.  However,  the 
governing  equations  for  a stiff  electron  beam  are  two  coupled 
second  order  differential  equations,  and  the  analysis  of  these 
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equations  is  considerably  more  complicated  than  that  of  an  electron 
beam  with  uniform  axial  drift. 
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